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Abstract 

We re-examine and extend a recent analysis which showed that the p—co mixing 
amplitude cannot be unambiguously extracted from the pion electromagnetic form- 
factor in a model independent way. In particular, we focus on the argument that 
the extraction is sensitive to the presence of any intrinsic ujj — > rnr coupling. Our 
extended analysis confirms the original conclusion, with only minor, quantitative 
differences. The extracted mixing amplitude is shown to be sensitive to both the 
intrinsic coupling loj — > tttt and to the value assumed for the mass of the p°-meson. 
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1 Introduction 



The p—uj mixing amplitude is traditionally extracted from the pion electromagnetic (EM) 
form-factor, F n (q 2 ), as measured in e + e~ — > 7r + 7r~ 0,|2[]. The non-perturbative strong 
interaction effects that produce the significant enhancement in the interaction around 
y/q 2 ~ 750MeV have been successfully parametrised using the vector meson dominance 
(VMD) model. To say any more about this is to undertake a study in intermediate 
energy QCD || and is not our purpose - we accept VMD as given. Our interest lies 
in the appearance of the isoscalar uo meson in the isovector, p° meson dominated, pion 
form-factor. The traditional treatment neglects any intrinsic coupling of the u to the two 
pion final state (i.e., that not proceeding through a p° meson) and hence assumes that the 
p—uj mixing amplitude is purely real. An argument due to Renard and others shows 
that, within certain approximations, the imaginary part of the p—u mixing amplitude is 
cancelled by the intrinsic u decay contribution to the pion form-factor. This argument 
was critically questioned recently || and argued to be unjustified. Here we carefully 
re-examine and extend these arguments. We confirm the central conclusion, finding, in 
addition, a considerable sensitivity to the value assumed for the p mass. 

2 The isospin pure basis 

Analysis of p — u; mixing starts with the "isospin pure" fields, pi and uj. (For brevity 
we will henceforth write p in place of p°). These are by definition exact eigenstates 
of isospin. The photon to hadron interaction can be conveniently described using a 
formalism in which the renormalised vector meson propagator is a 2 x 2 diagonal matrix. 
For the pion EM form-factor, in the pure isospin limit, we have no p — uj mixing and no 
direct ujj — ► titt coupling and hence 

F = ~( f f ) I DlpP H 

where e = |e|, f 1 y I is the coupling of the vector meson to the photon and g Pl7T7T is 
the coupling of the p to the two pion final state. Here D 1 ^ and are the scalar 
parts of the renormalised propagators for the isospin pure fields. As is standard in 
traditional treatments, coupling to conserved currents is assumed, allowing us to use 
D[Mv(q 2 ) = —g pv D(q 2 ), where D(q 2 ) is the "scalar propagator." It is important to note 
that unless one assumes that the vector mesons do couple to conserved currents, the q^q v 
pieces of their propagators must also be included in the analysis. 

To model the G-parity violation seen in F n (q 2 ) one needs to somehow include the 
uj pole in the n + ir~ production by breaking the SU(2) isospin symmetry of the system. 
The usual mechanism has the two eigenstates "mix" through the isospin violating mixing 
self-energy, Il puJ (q 2 ), to allow the decay uj — > pi — > tttt through off-diagonal terms in the 
dressed, isospin- violating vector meson propagator matrix 0. This introduces isospin 
violation (IV) in the vector meson propagator and we have (retaining only first order in 
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isospin violation) || 



D 1 = [ D w ° ] - D 1 




d 1 n iVid 7 

J ^pp 11 pLj\H I^LUL 

^ pp LL pw\H ) lulu 



+ 0((IV) 2 ).(2) 



LULU 



However, a priori, in an effective Lagrangian model involving the vector mesons, all 
isospin violation sources are equally likely and there is no reason to exclude the possibility 
of the "intrinsic decay", uj — > nn, through the coupling g Wl7rn . The appropriate VMD- 
based expression for the pion form-factor in the isospin pure basis would then be given 
by 

1 , . .J dL dIjl„,dL 



pp - L ^pp J - J -pU J 
6 \ ^ppLLpLU-^uJuj ^LULU 




fjpi 1 , f~yui 1 o / 2\ ^ 

2 2/ 2\9pI 7Tn ~\ 2 2/ 2\ V^W J - 2 2/ 2\9pii™ 

e q z — rrip{q z ) e q z — m^{q z ) q z — mp(q z ) 



2 _ rr,2( n 2\^I^' 



e q 2 -ml(q 



(3) 



where a fourth term on the RHS involving both and has been neglected since 
it is second order in isospin violation. One is always free to consider models where g WJ7r7r 
is strictly zero, but there is no model-independent requirement that this be so. For 
the renormalised, isospin-pure propagators, D vv , we have used the physical p and uj 
propagators, since D vv = D vv + C((IV) 2 ) || and since again we are working only to 
first order in isospin violation. For the physical p and u propagators we will use here the 
usual Breit Wigner form with a momentum dependent width, i.e., 

Dyv = ~2 rr^T' ( 4 ) 

q z — m v {q z ) 

where 

m v(l 2 ) = ™v ~ irh v T v h v (q 2 ) (5) 

and where we have defined the momentum- dependent width Ty(q 2 ) = Tyhy(q 2 ) with 
hy(rh v ) = 1 and where the mass parameter (my) and width parameter (IV) are real. 



3 The physical basis 

The pion form-factor has two distinct resonance poles, yet Eq. (§) has a rather compli- 
cated pole structure. Therefore, one might change to a basis that has only a sum of simple 
poles. This can be done by transforming to what is referred to as the "physical" basis. 
It has been traditional to assume that ILj^ is a constant which allows the off-diagonal 
terms in the propagator to be removed by a rotation of the fields. However, recent work 
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has argued that, when p—oo mixing is generated only through vector mesons coupling to 
conserved currents, this amplitude is necessarily momentum dependent 0. This result 
is consistent with an earlier study which concluded that momentum independent mixing 
was acceptable for scalar particles, but not for vectors coupling to conserved currents |7[. 
Various demonstrative model calculations support this ||. To make allowance for this, 
Maltman, O'Connell and Williams (MOW) defined the physical basis through two isospin 
violating parameters f5|, 

p= Pi- exLUj, 
uj = uji + e 2 pi- 

This defines a transformation matrix, C, via 






The off-diagonal element of the physical propagator was then defined in terms of the 
physical fields in the usual way through 

D%{q 2 ) = i J d 4 xe^(0\T(p^x)u;%0))\0). (8) 

As the vector mesons couple to conserved currents, we can replace D^ u {q 2 ) by —g^ y D{q 2 ). 
Eq. (Bf) leads us to 



DM) = D^(q 2 ) - e l DUl 2 ) + e 2 D I pp (q 2 ) 

= DjvtfWM) - ^(g 2 ))- 1 + e 2 (DUq 2 ))- 1 ]DL(<l 2 )- (9) 

Note that 

D = CD 1 ^. (10) 

The physical basis is defined by requiring that there are no resonance poles in D puJ (q 2 ), 
i.e., we choose e\ and €2 such that D puJ (q 2 ) has no pole. The possible pole positions are 
the vector meson pole positions, m 2 and m^. These singularities can be removed if the 
numerator of D puJ (q 2 ) vanishes at these positions, leading us to, 

— m p — m p 

It should be noted that m 2 and m^, are the complex resonance pole positions, as given 
in Eq. ([5]). In the case where the mixing is momentum dependent we have e\ 7^ e 2 such 
that the matrix C is not orthogonal and the transformation between bases is then not a 
simple rotation, i.e., 

cTc = ( £2 - £l V)*'- < 12 > 

In a similar fashion to Eq. (|8|), the coupling constants in the physical basis are defined 
via 

f^U) 4" tzfjpjl QljJIT-K QuJ^TT 4~ ^29piTT1T (13) 

fjp fjpl ^lf-yuljl 9p1TTT 9plTTTT- 
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Note that eig W/7r7r is second order in isospin violation and so is not retained. 
In this way, MOW defined the form-factor in the physical basis by 

F n(q 2 ) = ~ [flWirAtfu>/ 7 w + gpnirDppf-yp + Qfm^D puf^} + background. (14) 

The non-resonant term, D pu) , was then included in the background, along with any other 
non-pole terms of the Laurent series expansions of the propagators, to arrive at an ex- 
pression for the form-factor, 

F^q 2 ) = H(e 1 , e 2 ) [P p + A(e 1 , e 2 )e i ' t,{ ^ ) P w ] + background. (15) 

Here P PyLU are simple poles and H is an overall constant. The quantities H, A and the 
Orsay phase, 0, can be read off from Eq. fll4]) as we will show later. In terms of the 
isospin-pure basis and the transformation matrix, C, Eq. ( |i~4| ) can be written, 



F-k — ~{fipi fjw I )C T CD I C T C ^ 




(16) 



With Eq. ( |i~2|) in mind, the two expressions for the form-factor (Eqs. (|T4| ) and (|IED) 
are not identical when e\ ^ e 2 . However, because of the closeness of the p and u pole 
positions, MOW find there is no practical distinction in their analysis between the two 
terms [i.e., U pLU (m 2 ) and n pLJ (m^)], so we define a single parameter, 



e 2 = — P — ^ ei (17) 

m A , — mi 



where 



n pw (mj) ~ U P U m l) - u pu(™ 2 p ) - n P ^(m^ 



In other words, from this point on we will assume (as is done in all standard treatments) 
that the momentum dependence of H. puJ {q 2 ) is negligible in the vector meson resonance 
region. In more careful language, we simply absorb any momentum dependence of 11^ 
into the non-resonant background. The transformation matrix, C, is now orthogonal, and 
the two bases are strictly equivalent. Note that this is entirely consistent with Eq. (|T0|) . 
In the limit that e\ = e 2 we find C T = C -1 and we recover Eq. (|16|). 

As a brief aside, we note for completeness that if ei and e 2 were quite different, the 
two bases could still be related in an appropriate way. The generalisation is simply to 
start with Eq. (^), but introduce the transformation matrix, C, in the following way, 

F-k — ~ [f-rVi] [D 1 ] [gvmn] 

= i[/ 7y7 ][^][^]- 1 [ J D / ][C]- 1 [C][^] 

= \[f lV ][D][gv™l (19) 

where C is chosen to remove resonance behaviour in the off-diagonal elements of D and 
where D = D when e\ = e 2 . 
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4 The Renard Argument 



The major conclusion of the MOW analysis || concerns the competition between the 
two sources of isospin violation, namely, 11^ and g Wl1tn . All isospin violation in the pion 
form-factor is usually attributed to p— uo mixing 0, because the intrinsic decay of the u> 
is assumed to be cancelled [Q. MOW paid careful attention to this Renard cancellation, 
which we shall now briefly review. Consider the uo pole term of Eq. ([TJ]). The coupling 
of the physical u to the two pion final state is given by, 

+ tg Pl 7m, (20) 



where we are using the e of Eq. fllTI). It is certainly a reasonable approximation to assume 
that the two pion intermediate state saturates the imaginary part of U puJ (q 2 ) around the 
p resonance region and that this is proportional to the two pion piece of the p self energy 
(which dominates the imaginary piece of the total p self energy, due to the strong p — > ixix 
decay). These are also standard assumptions. We have then 

Imn^(mJ) = kn™(m;) 

= f^imn-K) 

ilpjITIT 

= -Gm p T p h p {m 2 ) 

~ —GrhpTp (21) 



where we have defined 



G = (22) 



9 P 



J7T7T 



and assumed h p (m p ) ~ h p (rh p ) = 1. We therefore define 

U pul (m 2 p ) = n pw (m 2 p ) - iGrh p T p . (23) 

Assuming saturation of the absorptive part by the two pion state, as described above, 
implies that U puJ (m 2 ) is real. Note that while the three pion state is kinematically acces- 
sible it also requires isospin violation and is, in addition, suppressed by the smaller phase 
space available. We follow standard practice and ignore it here. 
Substituting Eq. ( |23|) into Eq. fllTD , we have 



n pw (m2) ^ rh p T p 



m 2 , — m 2 m 2 , — m\ 



p 

= -izf-zG, (24) 



where we have defined 



irhpTp ~ _ U puJ (m p 



m 2 -m 2 p rhpYp 



2 ) 

(25) 
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Note that since in the resonance region we are neglecting momentum dependence of the 
resonance widths, i.e. we take h^q 2 ) = h p (q 2 ) = 1, then from Eq. (||) we see that m^, 
and z are all constants. Recall that here, since we are using just the single parameter 



ni 



e, the mixing amplitude H puJ (m 2 ) should be understood to mean H puJ (q ) with q in the 
vector resonance region. Using this expression for e in the uj pole term of Eq. (^) gives 
rise to Renard's cancellation 



JUJJTYTT 



£9 P 



J7T7T 



jLOjTTTT 



'l-z) 



G(l-z) + 



m 2 , — in 



^(m 2 p ) 



mi 



2 9pinin 



9pv 



(26) 



when one makes the approximation z — 1. However, MOW find that z has a sizeable 
imaginary piece 

z = 0.9324 + 0.3511 i, (27) 

using the mass and width values from Bernicha et al. |T(|. For comparison, one finds 
z = 1.023 + 0.2038i using the values of Benayoun et al. [fTI |. The central observation 
of MOW was to point out that the deviation of z from unity leads to a substantial 
contribution to F T (q 2 ) from g W/7r7r . 

This brief discussion summarises the essential arguments of MOW. At this point 
MOW neglected the e dependence of the constant, H(ei, e 2 ) in Eq. (|TS|) and extracted G 
and Yl puJ (m 2 ) by fitting A and the Orsay phase <fi. In order to do a more careful analysis 
of the G dependence of the extracted real part of the p — u mixing amplitude, fl pu j(m 2 ) 



we shall retain H(e) = H(ei, e 2 ), which itself has a G dependence [see Eq. (p2|)j through 
e [via Eq. (E3I)]. Let us write Eqs. (fTJ) and flT5|), [or equivalently Eq. (^)] as 



Fir(q ) ^[f / ypi9piTTirPp f^u)i^i,Pp Pw)9 P inn fjoji Puj9(jJitttt\ 



where we have used the fact that e(P p 
ratio 



P,, 



PpUpuPu using Eq. (|1 



77 = 



IPi 



(28) 

Defining the 
(29) 



we can use Eqs. 

FM 2 ) = 

Comparing Eq. 

identify the qualities H, A and the Orsay phase, <fi. This is the central result. 



H) and (H) to write Eq. (H) as 

-f 1Pl 9 Pl ™ {P P + r/ [-e(P p - PJ) + GPJ} 
-f 1Pl 9 Pl ™ {Pp{^ + izrif) + rj [-zzf P w + zGP p + G(l - z)P u 



} (30) 



to our earlier Eq. flT5|) it is seen to be a straightforward matter to 
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It is instructive to now collect together all of the terms inside the braces {• • •} in 
Eq. (p0[) containing G. We have 



Gr I [P u} (q 2 ){l-z)+P p {q 2 )z] 



G ri P w {q 2 )P p {q 2 ) 
G ri P w {q 2 )P p {q 2 ) 



1 



1 



\P P (q 2 ) W). 
z (m u - rrip) 



q 



m. 



(31) 
(32) 



in the resonance region. We see that in this complete analysis, the total contribution 



associated with G vanishes identically at the point q 



in 



regardless of the numeri- 



cal value of z. Recall that to obtain this result, we followed standard treatments and 
neglected the momentum- dependence of the Breit-Wigner widths in the resonance region, 



h p {m p ) ~ h p (rfi p ) = 1. 



(33) 
(34) 



We also neglected the three pion loop contribution to p — u mixing, so that 

ImlL^g 2 ) = GlmU pp (q 2 ) 

and assumed a constant value for Tl pL0 in the vector meson resonance region. 

We arrive at a similar conclusion in a more straightforward way if we choose to stay 
in the isospin pure basis throughout and simply re-arrange Eq. ([3]) to give 



fAq 2 



7P V 



m 2 (q 2 ] 



+ 



7Wj 



q 2 — m 2 (q 2 



■n p (g 2 



9pi 



q 2 -ml(q 2 



+ 



70)/ 



im p Y p h(q 2 



q 2 — m 2 J (q 2 ) 



JUJJ7V1T 



q 2 — rh 2 + irh p T p h(q 2 ) 



(35) 



We see immediately that in the isospin pure treatment the G dependence (i.e., g, 
dependence) is cancelled at q 2 



LUjTYTV 



m. 



and somewhat suppressed around the pole region 
(where p — u interference is most noticeable). Note that in writing Eq. fl3"5|) we have 
again made the assumption of Eq. (0). Finally, we note that Eq. (p2|) has been derived 
by ignoring the non-resonant contribution, D pLU (q 2 ), to the form-factor. Including this 
would make Eqs. (|30D and (|35|) identical, but as we can only fit resonant terms plus an 
unknown non-resonant background to data, there is no practical difference between these 
two equations. 

Summary: We have shown that when the standard approximations (i.e., neglecting 
momentum dependence of both the mixing self-energy and the Breit-Wigner widths in 
the resonance region) are imposed, then the results of MOW and those of the present 



work are exactly equivalent at q 



m. 



and should be very close near this point. It now 



remains to extend the numerical analysis of MOW by keeping all of the G dependence 
[i.e., including the G-dependence of H in Eq. flip] ) as given by Eq. (|30|)1. Assuming that 
the pole positions and isospin conserving couplings are already known, we are left with 
just two parameters, ^ J7r7r and ti^m 2 ) (or equivalently G and T) to be fitted to the 
data. 
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5 Results and conclusions 



To determine values for n^m 2 ) and G we must choose an appropriate form factor and 
decide upon Eq. fl35|) , which we re- write as 



—am. 



w? p + irhpTp 



1 + rj 



m. 



ml + im u T u 



(36) 



We now perform a fit to pion form-factor data to extract values for m p , T p , 



G 



and the normalisation constant, a, using the fitting routine PAW [|ll| (case A in Table [I]). 
We fix m u = 781.94 MeV and 1^ = 8.43 MeV, as given by the Particle Data Group |]|. 
It should be noted that we are here adopting this form to fit rather than using Eqs. (1ST 



and (15), since we wish to isolate all of the G-dependence from the overall normalisation 
constant for the fit. It can be seen from Eq. ( pCf ) that H in Eq. (|TJ) is G-dependent and 
so we can not use the fitting procedure adopted by MOW. By construction, the constant 
a in Eq. (^) is independent of G. For the analysis here we will assume the SU(3) value 
for n = g^ux/g-fpu i.e., we will assume 77 = 1/3 Hl3| . However, from Eq. (|36|) we see 
that all that can really be extracted from the analysis is riU pLU and rjG and so if another 
value of rj is preferred our results can immediately be scaled in a straightforward way. 
Our data set consists of the 70 points listed in Ref. fl4| in the region between 500 and 
975 MeV. While the preferred value for G is quite large - one usually expects isospin 
breaking at the few percent level, rather than 10% - the uncertainty is such that it lies 
only 2 1 standard deviations from zero. 

We can investigate the importance of direct isospin violation at the uo — > tttt vertex 
by imposing the condition G = (case B). However, when fitting the p data, it is 
important to remember that the p parameters should not be process-dependent. As noted 
by Benayoun et al. since the p is a relatively broad resonance the value extracted for, say, 
the mass, can be greatly affected by the addition of other terms to a phenomenological 



form-factor pT| . To demonstrate that our general conclusions concerning G and n pw do 
not depend on a particular choice of mass and width from one source, we have redone 
the fit using the Particle Data Group (PDG) values, which are averaged from a variety 



of processes [15 



In order to compare our results with the analysis of Maltman et al. [||], we perform a 
fit to the same data using Eq. (|TJ) with H(e) treated as a constant (thereby absorbing 
it into the normalisation constant a). To first order in isospin violation one has 



r(G(l - z)-izT) 
1 + izrT + zrG 



(37) 



The resulting fit parameters are shown in Table 1 under the heading MOW. We see that 
they are very close to the fit using the complete expression Eq. ([30]), as shown in column 
A. (Note that the errors on n and G are closely correlated, but this is reflected to a 
certain extent in the diagonal errors which are determined in the standard way through 
inversion of the curvature matrix.) 
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Parameter 


A 


B 


PDG 


MOW 


Il pu] (m 2 p ) (MeV 2 ) 
G 
a 

m p (MeV) 
T p (MeV) 
X 2 /dof 


-6832 ± 1252 
0.102 ±0.04 
1.15 ±0.02 
764.1 ±0.7 
145.0 ± 1.7 
88/65 


-3844 ± 271 
(input) 
1.19 ±0.01 
763.7 ±0.7 
146.9 ± 1.4 
94/66 


-6298 ± 788 
0.065 ±0.033 
1.203 ±0.004 
769.1 (input) 
151.0 (input) 
172/67 


-6827 ± 1067 
0.102 ±0.04 
1.20 ±0.01 
764.1 ±0.7 
145.0 ± 1.6 
88/65 



Table 1: Results from fitting our form-factor to data. 



In conclusion, we have seen that, in agreement with the conclusions of Maltman et al, 
the pion form factor data supports equally well a large range of possible pairs of values 
for G and Tl puJ . In other words, it is not possible to extract the p—uj mixing amplitude in 
a model-independent way. The traditional method of extraction corresponds to assuming 
that there is no intrinsic u — > titt coupling (i.e., that G — 0), which is highly unlikely. 
It should also be noted that these conclusions are entirely independent of what (if any) 
momentum-dependence is present in the p— u mixing amplitude, since this was neglected 
in the resonance region in the usual way. 
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